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Abstract 

In quintessential inflationary model, the same master field that drives in- 
flation becomes, later on, the dynamical source of the (present) accelerated 
expansion. Quintessential inflationary models require a curvature scale at the 
end of inflation around 10 -6 Mp in order to explain the large scale fluctuations 
observed in the microwave sky. If the curvature scale at the end of inflation is 
much smaller than 10 _6 Mp, the large scale adiabatic mode may be produced 
thanks to the relaxation of a scalar degree of freedom, which will be generically 
denoted, according to the recent terminology, as the curvaton field. The pro- 
duction of the adiabatic mode is analysed in detail in the case of the minimal 
quintessential inflationary model originally proposed by Peebles and Vilenkin. 
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I. INTRODUCTION 



Following the discovery [1,2] that distant supernovae are fainter than inferred from local 
samples, models of scalar fields that are able to develop a negative pressure at the present 
time have been proposed. These scalar fields only interact gravitationally and they have 
been generically named quintessence [3-7]. 

In recent years, an interesting class of quintessence models have been proposed by Peebles 
and Vilenkin [8] (see also [9]). The peculiar feature of quintessential inflation is that the 
master field that drove inflation becomes again dominant, at the present time. Hence, in 
the context of quintessential inflation, the inflaton and the quintessence field are identified 
in a single scalar degree of freedom, driving inflation in the past and acting as quintessence 
today. A simple example of this dynamics is represented by "dual" potentials [8] going as a 
power during inflation and as an inverse power after [3] inflation. Owing to this dual form of 
the potential, the dynamics of the background will experience, right after inflation, a pretty 
long phase in which the potential term is subleading with respect to the kinetic term. In 
quintessential inflation the reheating can be entirely gravitational, since the energy density 
in the light quanta produced at the end of inflation decreases slower than the (kinetic-energy- 
dominated) background geometry [8,10]. 

Quintessential inflationary models are consistent with observations provided the curva- 
ture scale at the end of inflation is not too low and O ~ (lCT 6 Mp). In the present paper a 
complementary possibility will be analysed. Consider the situation where H e , the curvature 
scale at the end of inflation, is indeed much smaller than 1CT 6 M P . Along this line we sup- 
pose, for simplicity, that during inflation there is degree of freedom, ip, which is not coupled 
with the inflaton field and remains constant, thanks to its potential, during the later stages 
of inflation. The (potential) energy density of ip is subleading with resepect to the energy 
density of ip. At the end of inflation, during the kinetic phase, the large scale fluctuations 
of ip will be converted into adiabatic fluctuations. This possibility has recently been studied 
by many authors [11-17] in different contexts, and it was originally invoked in [18]. As dis- 
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cussed in [19], even the simplest chaotic inflationary models develop new constraints when 
combined with the curvaton idea. The purpose of the present paper is to analyse low-scale 
quintessential inflation in a specific set-up, which is the one originally suggested by Peebles 
and Vilenkin in [8]. In this model the late-time behaviour of the quintessential evolution 
does not show the tracking behaviour of the inflaton and matter mass densities, as argued 
in [6]. 

Particular attention will be given to the evolution of the fluctuations. This is quite 
essential since, by lowering the curvature scale at which inflation ends, we have to make 
sure that the isocurvature mode of ip is efficiently turned into an adiabatic one. While in 
ordinary inflationary models (such as those analysed, for instance, in [19]) the radiation- 
dominated phase starts at the end of inflation, in quintessential inflation the evolution may 
be very different and the onset of the radiation-dominated phase is delayed. One of the 
purposes of the present paper is indeed to generalize the analysis of the curvaton evolution 
to backgrounds where inflation is not immediately followed by a radiation-dominated phase. 
The curvature perturbation will be followed through all the stages of the model and its final 
value computed. This analysis will be performed both analytically and numerically. 

The plan of the present paper is the following. In Section II the constraints on the post- 
inflationary evolution will be derived in the specific case where the inflaton and quintessence 
field are identified. In Section III the constraints pertaining to the quintessential evolution 
will be scrutinized. Section IV contains the basic ingredients for the evolution of the fluc- 
tuations. Sections V and VI deal with the conversion of the initial isocurvature mode into 
an adiabatic one. In Section V the initial and the kinetic stages will be scrutinized, while 
Section VI is more oriented towards the phase where ip dominates and eventually decays. 
In Section VII some concluding remarks will be presented. 
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II. FROM INFLATION TO QUINTESSENCE 



M 2 H 2 



* +t- + V(<p) + W(rl>) , (2.1) 



Consider the minimal realization of a low-scale quintessential inflationary model extend- 
ing the original proposal of [8] : 

L2- + T 

M 2 (H 2 + H)= [-y? 2 -ip 2 + V((p) + W(ip)\ , (2.2) 
dV 

ip + 3Hip+ — = 0, (2.3) 
dip 

dW 

ip + 3Hip + — - = 0. (2.4) 

dip 

The potential of ip can be chosen to be a typical power law during inflation and an inverse 
power during the quintessential regime: 

V(<p) = A(v? 4 + M 4 ), V <0, 
A /If 8 

where A is the inflaton self-coupling and M is the typical scale of quintessential evolution. 

The field ip is subleading during inflation and it is characterized by a potential, which 
we will take, for simplicity, to be quadratic, i.e. 

o 

W(ip) = —ip 2 . (2.6) 

This set-up can be generalized to the case where the field ip is replaced by an arbitrary 
number of scalar degrees of freedom ipi. 

Inflation ends when V((p) ~ AMp at a curvature scale H e 

In ordinary inflationary models, the large scale inhomogeneities determining the CMB 
anisotropics come from the Gaussian fluctuations of the inflaton ip, and, consequently, 
A ~ 1CT 13 [8]. In the present investigation a complementary possibility will be discussed, 
namely the case when 



H c < KT 6 M P . 



(2.8) 



In this case the fluctuations of <p are too small to be interesting for CMB physics. However, 
this conclusion can be evaded by taking into account the fluctuations of the field ip. Quali- 
tatively the picture is the following. Right after inflation, the fluctuations of the geometry 
will be determined by the fluctuations both of (p and ip. Because of Eq. (2.8), the metric 
fluctuations generated by ip will vanish at t e . In other words, the initial conditions for the 
system at t e will be of isocurvature type. Later on, the isocurvature mode can be converted 
into an adiabatic one. This is the idea explored, for instance, in [11,12]. In the usual picture 
of curvaton evolution, right after inflation, radiation takes place immediately. In the case 
of quintessential inflation, on the contrary, the onset of the radiation-dominated epoch may 
be delayed. Hence, the analysis of the evolution of the fluctuations must be repeated, with 
particular attention to the different features of the model. 

During inflation, the field ip is subdominant with respect to the inflaton energy density, 



The field ip remains nearly constant during the later stages of inflation, i.e. ip ~ ip e . At the 
end of inflation, Eq. (2.9) implies that 



After the end of inflation, because of the inverse power-law form of the potential, the 
field ip is mainly driven by its kinetic energy and the approximate solution of the background 
geometry will be of the type 



with a(t) ~ t 1 / 3 . In this phase the potential, V(ip) ~ A(M/M P ) 4 In -4 (a/o e ), is subleading 
since A <C 1CT 13 and M <C M P . More specifically, phenomenological considerations related 
to the present dominance of (see Section III) suggest that M <C lCT 9 Mp. 



W(iJ>) < V((p). 



(2.9) 




(2.10) 




(2.11) 
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FIG. 1. The result of the numerical integration of the background are illustrated for different 
values of the parameters (indicated above each curve) . 

While the background is dominated by <fi 2 , the field ip slightly decreases with a rate given 
by W^. From Eq. (2.4), the field if) slowly rolls toward the minimum of its potential in a 
kinetic-energy-dominated environment and its approximate equation obeys 

1 dW 



(2.12) 



QH dip ' 

where the factor 1/6 comes from dropping consistently the terms of Eq. (2.4) containing 
more than one derivative of the potential. This type of approximation has been exploited 
in [13], but in the case of a slow-roll occurring during the radiation epoch. 

Using Eq. (2.12), and recalling that a(t) ~ i 1 / 3 , it can be checked that the evolution of 
ip slightly deviates from a constant value. For instance, in the case of a quadratic potential 
the solution of Eq. (2.12) can be written as 



m , 9 

1 (t 2 

4 v 



(2.13) 



Equation (2.12) is useful in the case when the potential is more complicated than the 
quadratic ansatz of Eq. (2.6) where Eq. (2.4) can be solved exactly. In fact, during 
the kinetic regime the curvature scale decreases, and when H ~ H osc ~ m the field ip, still 
subdominant, will start oscillating. The case of massive potential allows analytical solutions 
of Eq. (2.4) during the kinetic phase for ip: 



* = *7*&L) J °W- (2 ' 14) 

which leads to Eq. (2.12) in the limit mt <C 1 and mt c <C I. Since ip decays as a~ 3 / 2 during 
the oscillating regime, it will become dominant with respect to the kinetic energy of ip at a 
scale 



The field ip eventually decays at a curvature scale 



H, ~ J. (2.16) 



At if m the potential of ip is still subdominant with respect to the kinetic energy. In fact 

^,^(^^(|)<f4(M )> ( , 17) 



where the last equality comes from the kinetic term of ip at t m , evaluated on the basis of 
Eq. (2.11). An example of the numerical integration is reported in Fig. 1. Furthermore, 
the analytical evolution of ip, as reported in Eq. (2.14), is in excellent agreement with the 
numerical results. 

Between H<± and H r , 

(?) - (*)'• w 

the background geometry is effectively dominated by the oscillations of ip, i.e. a(t) ~ (mt) 2 / 3 . 

It is interesting to combine in the physical picture the constraints and the requirements 
introduced so far. According to Eqs. (2.15) and (2.16) the quintessence field may become 
subdominant either prior to or after the decay of ip. In Fig. 2, with the (diagonal) dashed 
line, the condition coming from the interplay between the decay of ip and the dominance of 
ip is illustrated for the specific case H e ~ 1CT 9 M P . Above the dashed line, ip decays during 
the kinetic phase. Below, the dashed line the decay occurs when already dominates. The 
requirement of Eq. (2.10) also imposes a constraint on Fig. 2, implying that %p e and the 

6 




-18 -16 -14 -12 -10 -8 -6 -4 -2 

log(V e /M p ) 

FIG. 2. The quantitative constraints pertaining to the combined analysis of the evolution of ip 
and ip during the kinetic phase are illustrated. With the dashed line the constraint coming from 
the dominance of ip is reported. The dot-dashed lined (from top to bottom) represent, respectively, 
the bound (2.20) and the BBN bound. The full (vertical) line refers to Eq. (2.19). The shaded 
region defines the allowed portion of the parameter space of the model. 
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mass should lie below the full fine in the right-hand corner. Consider now the situation 
where the fluctuations of if) are amplified, during inflation, with a scale-invariant spectrum. 
Therefore, denoting by Xi> the fluctuation of if), we do know that the power spectrum will 
be 5 M ~ H c /2n. If 

TT < 27l TT> ( 2 - 19 ) 
M P M P v ' 

then the fluctuations of if) will be predominantly Gaussian. In the opposite case they will 
have some non- negligible non-Gaussian component. In order to have a nearly scale-invariant 
spectrum for Xif>, m should always be smaller than the curvature scale H e . More precisely, 
looking at the evolution equation for Xtp (in conformal time) we are led to require 

m < V2H C . (2.20) 

Equation (2.20) is illustrated in Fig.2 with the dotted horizontal line. Finally, since the 
decay of if) should occur prior to BBN, an absolute lower bound on m, i.e. m > 10 TeV, 
should be imposed. Therefore, already from these considerations it is possible to say that the 
dashed area in Fig. 2 is allowed. Thus, the quintessence field has to become subdominant 
before the decay of if) if we want Xtp to be Gaussian with nearly scale-invariant spectrum. 



III. QUINTESSENTIAL EVOLUTION 

In order to satisfy the physical constraints for the consistency of the scenario, it is 
reasonable to require that the decay of if) occur when if is already sub dominant. If this is 
the situation, > H T and, consequently, 

lf) e ~"~^ m - ( 3 - 1 ) 

It is important, for the purposes of the present investigation, to analyse in some detail, the 
evolution of ip around the curvature scales and H r . 

After Hd, the evolution of tp results from the interplay between the smallness of its 
potential term V(ip) and the coherent oscillations of if;, which tend to make the average 
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expansion matter-dominated. The effect of the potential is, however, to introduce a growing 
mode in the evolution of ip. This growing mode can be estimated by solving Eq. (2.3) with 
the inverse power-law potential of Eq. (2.5). The solution is 

p g (t) ~ 3 1 / 3 A 1/6 M 4//3 t 1/3 , t d < t < t T . (3.2) 

For t > t r , when the Universe is effectively radiation-dominated, a similar solution holds: 

^(t) ~ (^) A l/6 M 4/3 f l/ 3> t > ^ (3 3) 

The effective evolution of p depends upon the balance of the growing mode with the 
decaying modes obtained from Eq. (2.3) when the potential is neglected. During the regime 
where tp dominates and in the subsequent radiation-dominated regime, the evolution of p is 
given, respectively, by 



W (t)^{[ln(|)+l]-|}, t,<t<t„ (3.4) 

Equations (3.4) and(3.5) are continuous in t T . Furthermore Eqs. (2.11) and (3.4) are con- 
tinuous in t d . As obtained in [8] the growing solution never comes to dominate against the 
decreasing mode. If 

^<)(^) 4 >(ir- < 3 - 6 » 

the growing mode becomes dominant only after the present expansion time t ~ Hq 1 . We 
will see that this condition is always verified because of the smallness of A. 

After t r the quintessence field is nearly constant and its potential energy is given by 

The condition that the potential energy in the quintessence field is comparable with the 
present energy density, 

V{f<) = Hf s Ml, (3.8) 



v ^ x mHz)V- (3 ' 7) 



fixes a relation among M and A 



M P 



M 




(3.9) 



Note that, for the typical parameters discussed so far, M is always greater than 10 5 GeV. 



Having discussed the main features of the post- inflationary evolution, it is now mandatory 
to understand the behaviour of the fluctuations. Available calculations on the curvaton 
dynamics in different models [11-17] always deal with a post-inflationary phase, which is 
dominated by radiation. Here, as discussed in the previous sections, the scenario is different. 
Right after inflation the background evolution is not determined by radiation but by the 
dynamics of ip, whose potential is now very steep. The fluctuation of the metric induced 
by the fluctuations of ip will be, at the onset of the post-inflationary epoch, very small 
since H e < 10 _6 M P . Hence, had we to classify the metric fluctuations at t C) we would say 
that their modes are of the isocurvature type. However, as time goes by, the curvature 
fluctuations, initially negligible at t e , will be driven to a constant value so that, much later, 
the initial isocurvature mode turns into an adiabatic one. 

In order to describe this picture quantitatively, we are led to consider the coupled system 
formed by the fluctuations of the metric and by the fluctuations of p and ip . The fluctuations 
will then be discussed in the longitudinal gauge, where it is particularly simple to relate the 
gauge-dependent quantities to gauge- invariant observables [20,21]. 

In the longitudinal gauge, the non-vanishing entries of the perturbed metric are 



where we used the fact that the fluctuations of the energy-momentum tensor are free of 
shear. For the fluctuations of p and ip, in the longitudinal gauge the following notation will 
also be adopted: 



IV. LARGE SCALE FLUCTUATIONS 



Sg 00 = 20, Sgtj = -2a 2 (p, 



(4.1) 
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y?^y9 + <5y9, 8(p = x<p, 
With these notations the evolution equations of the fluctuations can be written as 



(4.2) 



+ 4#0 + (2H + 3H 2 )(j) 



2M 2 



dV 



3H(H<P + 0) 



a" 



,-V 2 = 



2M 2 



(if + if;) - (ip X ^ + VXv) + -q^Xv + 

dV 



dW 



dtp 



■Xi> 



dip 



dW 

dtp 



■XV 



i 9 <9 2 \/ • dv 

X„ + 3ff* v " ^V 2 X „ + -^X* - 400 + 2—0 = 0, 

1 d 2 W ■ ■ dW 

X<p + ZH X ^ - ^V 2 X ^ + ~^X^ ~ 4^0 + 2^0 = 0, 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



where (4.3) and (4.4), come, respectively, from the and (0,0) components of the per- 
turbed Einstein equations and (4.5) and (4.6) describe the evolution of the inhomogeneities 
in the inflaton/quintessence field and in the curvaton field. Equations (4.3)-(4.6) are sub- 
jected to the momentum constraint 

3 



#0 + 



2M 2 



(VX^+^Xi:)- 



(4.7) 



Using Eq. (4.7) together with Eqs. (4.3)-(4.6), it is possible to obtain a nicer form of the 
perturbation equations [22-24,26]: 



+ 3Hv v -V + 

a 



d 2 V 3_ d /a 3 .. 

~ Mjc 

d 2 W 



a fa J 9 \i 3 d /a 3 -\ 



v M 2 



3 d (a 



nTT . i ^ 2 \o*w 6 a tar ;2 \i 6 a (or -\ n , A >. 

v * + 3Hv " -* Vv * + lw~ m^diVH* )\ v " ~ m^diKH^h = °' (4 - 9) 



_3_ d /a 3 



where 



v< fi = X v + jj<P, 



(4.10) 
(4.11) 



Eqs. (4.8) and (4.9) can be generalized to the case of an arbitrary number of fields [22,23]: 
in this case the number of equations will clearly match the number of fields but the relative 
structure of the equations will be the same. 
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Eqs. (4.4)-(4.6) or, equivalently, Eqs. (4.8) and(4.9), have to be studied and solved 
along the different stages of the evolution of the background. In terms of v 9 and the 
spatial curvature perturbation can be written as 

C = ~ . 2 , ; 2 [foV + #^]- (4-12) 
The variable £ is related also to by the usual expression 

C = §(#0 + 0)-0. (4.13) 

Eq. (4.13) can be obtained from Eq. (4.12) (or vice versa) by using the momentum constraint 
(4.7) expressed in terms of the of v v and iy given in Eqs. (4.10) and (4.11). 

An interesting (complementary) strategy in order to solve the system (4.3)-(4.6) and 
(4.7) is to write down directly the evolution equation for (. Mutiplying Eq. (4.4) by the 
sound of speed and subtracting it from Eq. (4.3) we obtain, at large scales [25], 

d( H 



dt tp 2 + tp 2 

where 8p na A is, in our case, 



5Pnad, (4.14) 



<fynad = (C 2 g - 1)0(0 + V>) + (1 - (%)(<j>X V + fe) - (1 + Cl) + ^X^) , (4-15) 

with 

2 V , 2 fdV . 3W ;\ 

c^ = t = 1 + b^ + ^T0 • 4.16 

P 3#(y? 2 + ^ 2 ) W d ^ ' 

In the first equality of Eq. (4.16) p and p are, respectively, the total pressure and energy 

density of the system written in terms of the two background fields (p and -0. Notice that in 

order to get to (4.16) the background equations of motion (2.1)-(2.4) have been used. 

The initial conditions of the system (4.4)-(4.6) after the end of inflation are dictated by 

the smallness of A. Going to Fourier space and considering only the super-horizon scales, 

the initial conditions of the system are, for A <C 10~ 14 , 

0(Mc) = O, X <p(k,te) = 0, x4k,t c ) = ^- (4.17) 

12 



In terms of v v and v$, we have at the end of inflation, from Eqs. (4.8) and (4.9), on 
super-horizon scales, 

v v (k,t e ) = 0, tty(Me) = x>(Me), (4.18) 

as is made clear by inserting Eqs. (4.17) into Eqs. (4.8) and (4.9). Physically Eqs. (4.17) and 
(4.18) guarantee the absence of adiabatic modes at t e . This aspect can also be appreciated 
by looking at the initial conditions for ( 

C(Me)=0, (4.19) 

as they follow from Eq. (4.12). 

V. EVOLUTION DURING THE KINETIC PHASE 

Using the initial conditions given by Eqs. (4.17)-(4.19), the evolution of the fluctuations 
can be solved in slightly different but, ultimately, equivalent ways. In the following, as a first 
step, the asymptotes for the evolution of the fluctuations will be discussed analytically in the 
vicinity of t e . The consistency of the analytical solutions with the results of the numerical 
integration is an important check to be done. Following recent techniques developed in a 
different framework [13], the solutions in the vicinity of t e can be obtained without specifying 
the form of the potential for the field ip. Later on, in order to perform the numerical 
integration, the case of massive potential will be mainly discussed. 

A. The initial stages around t c for generic potential 

During this phase the field ip slowly rolls down its own potential W(ip), and the solution 
to Eq. (2.4) is given by (2.12). The same expansion in the gradients of the potential W(ip) 
can be used in order to get the approximate evolution of the fluctuations. 

Eqs. (4.8) and (4.9) can then be approximately solved by keeping the leading terms in the 
derivatives of W(ip). As done in the case of the background evolution, an equation analogous 

13 



to (2.12) can also be obtained for the canonical perturbation variable v^. Eq. (4.9) can then 
be expanded in gradients of the potential leading to the following approximate equation 

1 d 2 W 

whose solution implies that is approximately constant. For instance, inserting into Eq. 

(5.1) the explicit expression for W(ip) given in Eq. (2.6), we get 

V^(k,t) ~7ty(Me)[l- — (t 2 - ^e)J ? t<t m . (5.2) 

This expression has been verified numerically for various values of the physical parameters. 
Using Eq. (5.2), Eq. (4.8) can be solved in the same approximation. Note, in fact, that 
in Eq. (4.8) the term in square brackets is sub leading with respect to the other terms for 
two separate reasons. First V }lflf is negligible in its own right, given the smallness of A and 
recalling that M/M P is O(10 -13 ). Second, since a 3 (p 2 /H is constant, its time derivative 
appearing in Eq. (4.8) is also negligible. The approximate evolution of v v (k,t) is then 
given, at large scales, by 

v, + 3Hv, -^|(^)^ = 0, (5.3) 

whose solution, using Eqs. (2.12) and (5.1), becomes 

MM) * -4^(|)^(*»*)[o 6 (*) -Ate)}- (5.4) 

Given that during the kinetic phase tp ~ H, then v v (k,t) oc a 6 . Again, in the case of Eq. 
(2.6), from Eq. (5.4) we obtain 

v v (k, t) ~ (j^)m 2 (t* - tDx^k, Q, (5.5) 

where, according to Eq. (4.18), v 1 p(k,t e ) = Xip(k,t e ) has been used. With the results of Eqs. 

(5.2) and (5.5), from Eqs. (4.10)-(4.12) we can also obtain the approximate form of the 
evolution of (f>k(t), Ck(t) and x^(k,t): 
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3 dW 

Ck{t) ~ 2Ml^p X ^ Q[a{t) ~ a (5 - 6) 
Mt) - ~^Ck(t), (5.7) 

X,(k, t) - (|) ^xAAt) - a«(Q]. (5.8) 

The same results derived in Eqs. (5.6)-(5.8) and based on the solution of Eqs. (4.8) and 
(4.9), with the initial conditions dictated by Eq. (4.18), can be obtained by integrating Eqs. 
(4.4)-(4.6). As a check of the consistency of the approach, it is in fact useful to insert Eqs. 
(5.7) and (5.8) back into Eqs. (4.3)-(4.6) and see that they are satisfied during the kinetic 
phase and prior to the oscillations of ip. Furthermore, the evolution of (k, as obtained in 
Eq. (5.6), can be also obtained directly from Eq. (4.14). In fact, during the kinetic phase 
ip <C if and, from Eq. (4.16) evaluated in the kinetic limit: 

2 1 dW 

° s ~ 1 " 18M|# 4 dij ' (5 ' 9) 



Hence, from Eq. (4.15) we find 



dW 

5pnad = (5-10) 



Inserting (5.10) into Eq. (4.14) and performing the integral, we get exactly Eq. (5.6). The 
analytic expressions derived so far allow a full control of the initial conditions of the system 
in the vicinity of t e . 

B. Evolution for t c < t < id 

After the onset of the kinetic phase at t e , but before the dominance of ip at t d , the 
evolution of the system can be solved numerically for various sets of initial conditions; an 
example of this behaviour is reported in Fig. 3 for the case of the potential given in Eq. 
(2.6). 

It will now be shown that the numerical evolution can be very accurately reproduced 
analytically by direct integration of the evolution equation in a well defined approximation 
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m = 1CT J H e 
y = 0.01 M p 
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FIG. 3. The result of the numerical integration for the evolution of the fluctuations are illus- 
trated for the case M/Mp = 10~ 13 and for a set of fiducial parameters chosen within the shaded 
region of Fig. 2. In the left plot the analytical results (dashed lines) obtained in Eqs. (5.14), (5.17) 
and (5.19) are compared with the numerical ones (full lines). 

scheme. The evolution of the fluctuations will be analysed first using Eqs. (4.8) and (4.9) 
and then using directly the evolution equation for (, i.e. Eq. (4.14). 

In the case of the potential (2.6), Eqs. (4.8) and (4.9) can be written as 
d 



3V2 2 



Mi 



+ 3Hv^ + m 2 



1 + 



6 



3v/2 2 



o. 



(5.11) 
(5.12) 



where Eqs. (2.11) and (2.4) have been used. Recalling now the exact solution for the 
evolution of (2.4), i.e. Eq. (2.14), it can be easily checked that, in Eq. (5.12): 
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il) e m < 1 



(5.13) 



for t < ~ m 1 (Mp/?/' e ) 2 . Neglecting the term containing v v in Eq. (5.12), the solution 
for the evolution of is given by 

(5.14) 



J {x e) 



where x = mt and Jq{x) is the Bessel function of order [27]. Eq. (5.14) reproduces exactly 
the numerical solutions for reported, for a particular case, in Fig. 3. Since it is always 
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true that mt e <C 1 for the constraints displayed in Fig. 2, the interesting limits of Eq. (5.14) 
are for mt <C 1 and mt >> 1. In the limit mt <C 1 Eq. (5.14) reproduces, as it should, the 
time dependence obtained in Eq. (5.2). For mt > 1 we have [27], from Eq. (5.14): 

v^(k,t) = tty(Me)y^cos (x - ^ (5.15) 

Note the similarity between Eqs. (5.14) and (2.14), which is a simple consequence of the 
quadratic form of the potential. 

Inserting now Eq. (5.14) in Eq. (5.11), and integrating a first time between t e and a 
generic time t we get 

dv v _ 372 ^(fc,t c ) r / x F(x c ) i 

~^ ~ "2M P j ( Xc ) i x v° {x) + Ji {x) ) ~ ~r\ ' (5 - lb) 

where F(x e ) = x 2 c [J (x c ) 2 + Ji(x e ) 2 }. 

Direct integration of Eq. (5.16) implies that 

V v (k, X ) = _^g ^MMe) \ x 2(j2 {x)+2J 2 {x) _ J o{x)Mx) \- G{Xe) -2F(x e )l (5.17) 
4 J {Xe) L V J J 

where G(x e ) = x 2 [J^{x c ) + 2Jf(x e ) — J (x c )J 2 (x e )] and where, as usual, x e ~ m/H e . Taking 
the limit for large x of Eq. (5.17) the following result can be obtained 

v v (k,t) = -^^ip e v^(k,t e )(mt) + ^^ij e v^,(k,t c )(mt c ) 2 \n^, (5.18) 

7T 4 t e 

showing off the linear growth of v v with time. 

Inserting now the solutions for ty, and v v given, respectively, in Eqs. (5.14) and (5.17) 
into Eq. (4.12), the evolution of ( for t < t d , turns out to be, for mt > 1, 



Note that, for t = t d , Eq. (5.19) gives 

c(Mi) ^p|M). (5 , 0) 

Using the result of Eq. (5.19) into Eq. (4.13) the evolution for the large-scale modes of the 
metric fluctuations can be obtained: 
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Again we verified that the same results can be obtained directly by integrating the Hamil- 
tonian constraint (4.4). 

The time has come to compare the accuracy of the analytical expressions derived in Eqs. 
(5.14), (5.17) and (5.19). These expressions have been plotted in Fig. 3 (left plot) with 
the dashed lines. For the same set of parameters (and with the full line) the outcome of 
the numerical integration has been reported. In the right plot the results for a different 
set of parameters are displayed. The analytical expressions match rather accurately the 
numerical results. Notice that the small wiggles in the evolution of ((k,t), modulating the 
linear growth, are due to the fact that we decided to plot the full expression of ((k, t), which 
contains Bessel functions, and not only its asymptotic limit for mt > 1. The same evolution 
for ((k,t) derived in this section can be directly inferred from Eq. (4.14), recalling the 
approximate form of 5p na d- This calculation is reported in the appendix. 

It is now appropriate to compare the situation of low-scale quintessential inflation with 
the situation occurring in the more conventional case of curvaton models where the infla- 
tionary phase is suddenly followed by the radiation-dominated phase. In this case the ((k, t) 
variable grows as \fi (i.e. linearly in conformal time) before the curvaton becomes domi- 
nant. Here we found that the growth is linear in cosmic time. In spite of this difference, 
the final amplitude of (,{k,t) is given approximately by Xv>(^> ^dVVv m fact, in the case of 
a quadratic potential evolving in a radiation-dominated environment, 

where we now have a(t) ~ \fi. Integrating once the previous formula we get ((k,t) ~ \ft. 
If the evolution occurs during radiation, the curvaton will become dominant at a typical 
curvature scale H d ~ m^e/Mp) 4 [13]. Using this result, the amplitude of ((k,td) is given, 
as previously anticipated, by x^(k,td)/ip e - 
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VI. DOMINANCE OF V 



Since if; starts dominating the background at t d , for t > t d the evolution of the system is 
described by the following set of equations 



H 2 Ml 



^ , 2,2' 
— + m ip 



(6.1) 



EMl = (6.2) 
$ + 3Hip + m 2 ip = 0. (6.3) 

For large times Eqs. (6.1)-(6.3) lead to an effectively matter-dominated phase where the 
oscillations of tp, 

, . . , . . cos mt 

induce oscillations in the Hubble parameter and in the scale factor, which increases, on 
average, as 

It is not difficult to see that, under these conditions, the evolution of v$ is dominated by 
a constant mode. In Eq. (4.9) the term containing is always suppressed for large times 
since,from Eq. (3.4) 



^ ^M P |, (6.5) 



and the quintessence field goes very rapidly to a constant. Thus, in this regime, the evolution 
of v^p is given by 

v^(k,t) = v^(k,t d )—. (6.6) 
Eq. (4.8) allows us to deduce that 

v v (k,t)~r 3 . (6.7) 

Inserting now Eqs. (6.4)-(6.7) into (4.12) we find that, for t > t d , ( is frozen to its constant 
value. Therefore, right before the decay of ip, the metric fluctuation will be given by 
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0(M r )~^^. (6.8) 

After t r the field ip decays and the evolution of the generated adiabatic mode becomes 
standard, namely we have the fluctuations of the quintessence field evolving in a radiation- 
dominated environment together with the constant mode of (. The evolution equations for 
the fluctuations of the quintessence field will then obey, at large scales, 

d 2 V ■ dV 

X v + 3#x v + -q^Xv = 4 V^> - 2 ( 6 - 9 ) 



The other equations are the standard ones, namely 



3 „ 3 



2Mf 2Mp 



dV 



3H(H(f> + 0) = ^pA + —p> -00 2 + <PX V + jtX* , (6-10) 



dip 



+ 4H<P + (3H 2 + 2H)4> = ^pA - ^#0 ~X^ + (6.11) 

3 r A 

H<t> + 4> 



2M| 



^PX^ + ^PrUr], (6.12) 



5 r -40 = O, (6.13) 
w r - - = 0, (6.14) 

where <5 r = 5p r / p T and w r is the velocity potential. 

As discussed in the context of the quintessential evolution of the background, the 
quintessence field, in the present model, does not lead to a tracking behaviour, as also 
noticed in [8] . Thus, the evolution of the fluctuations during the radiation-dominated stage 
of expansion will effectively be the one implied by a standard cosmological term. In fact we 
can recall, from Eq. (3.5), that ip approaches a constant value as t -1 / 2 while the potential 
is constant. Then, from Eq. (6.9) it can be deduced that also ~ t~ l l 2 at large scales. As 
a consequence, combining Eqs. (6.10) and (6.11) we have, at large scales, 

+ 5#0 + 2(# + 2# 2 )0 = 0. (6.15) 

leading to the usual constant mode which was present prior to matter-radiation equality, a 
known feature of these types of models [28]. 
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VII. CONCLUDING REMARKS 



In the present paper, low-scale quintessential inflationary models have been analysed. 
The set-up of the model is closely related to the one proposed by Peebles and Vilenkin. 
Even if the curvature scale at the end of inflation is much smaller than 10~ 6 Mp, curvature 
perturbations of correct amplitude can be generated. In these models, the quintessence 
and the inflaton field are identified. The consistency of the background evolution imposes 
a number of constraints on the various parameters. Assuming a specific form of the in- 
flaton/quintessence potential and a massive curvaton, the evolution of the fluctuations has 
been followed through the different stages of the model. Even if, for some region of the 
parameter space of the model, non-Gaussian (adiabatic) fluctuations can be generated, the 
case of the Gaussian adiabatic mode has been discussed. The curvaton evolution occurring in 
low-scale quintessential inflation is different from the one possibly obtained in the case when 
the radiation-dominated evolution follows immediately inflation. In low-scale quintessential 
inflation, prior to curvaton dominance, curvature perturbations grow linearly in cosmic time. 
On the contrary, if inflation is immediately followed by a radiation-dominated stage, then 
curvature fluctuations grow as \ft. 

In spite of the fact that low-scale quintessential inflationary models share essential analo- 
gies with the prediction of ordinary quintessential inflation there are also relevant differences 
between them. In ordinary quintessential inflationary models, the reheating is mainly grav- 
itational. In the present case the reheating is triggered by the curvaton decay. In ordinary 
inflationary models a large background of gravitational waves should be expected in the 
GHz region [29,30]. In the case of low-scale quintessential inflation, gravitational waves are 
negligible over all the frequencies of the spectrum. 
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APPENDIX A: EVOLUTION OF ( DURING THE KINETIC PHASE 

In this appendix the analytical solution of Eq. (4.14) will be obtained for t e < t < t d . 
During the kinetic phase, ip 2 ^> ip 2 . Furthermore, as discussed in the paper 

dW 



Thus, from Eq. (4.14), we obtain that 

§ ~rj"'" (A.2) 

at (p z 



For t c < t < t d , the approximate evolution of ip and Xtf> is given by 

y/t Jo(mt) 



a 3 / 2 JqK)' 



X^)=X#,g^^. (A.3) 
Inserting these solutions into Eq. (A.2) we get 

c{k > t] = \ (jO (^t 1 ) M 2 ( J ° 2(mt) + J2(mt) ) - (mtc)2 ] • (A - 4) 

where we used the fact that mt e <C 1 because of the constraint stemming from Eq. (2.20). 
Taking the limit of Eq. (A. 4) for mt > 1 we obtain 



which is the same as Eq. (5.19) since at t c , according to Eq. (4.18), Xij>(k,t e ) = v^(k,t e ). 
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